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Abstract 

We have recently proposed an "Ultra-Strong" version of the Equivalence Principle (EP) 
that is not satisfied by standard semiclassical gravity. In the theory that we are conjecturing, 
the vacuum expectation value of the (bare) energy momentum tensor is exactly the same as 
in flat space: quartically divergent with the cut-off and with no spacetime dependent (sub- 
leading) terms. The presence of such terms seems in fact related to some known difficulties, 
such as the black hole information loss and the cosmological constant problem. Since the 
terms that we want to get rid of are subleading in the high-momentum expansion, we at- 
tempt to explore the conjectured theory by "IR-completing" GR. We consider a scalar field 
in a flat FRW Universe and isolate the first IR-correction to its Fourier modes operators that 
kills the quadratic (next to leading) time dependent divergence of the stress energy tensor 
VEV. Analogously to other modifications of field operators that have been proposed in the 
literature (typically in the UV), the present approach seems to suggest a breakdown (here, 
in the IR, at large distances) of the metric manifold description. We show that corrections 
to GR are in fact very tiny, become effective at distances comparable to the inverse cur- 
vature and do not contain any adjustable parameter. Finally, we derive some cosmological 
implications. By studying the consistency of the canonical commutation relations, we infer a 
correction to the distance between two comoving observers, which grows as the scale factor 
only when small compared to the Hubble length, but gets relevant corrections otherwise. 
The corrections to cosmological distance measures are also calculable and, for a spatially flat 
matter dominated Universe, go in the direction of an effective positive acceleration. 
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1 Introduction 



Modifications of General Relativity (GR) on the largest scales have been advocated in order to 
give account for the present acceleration of the Universe. Attempts in this direction [1] typically 
maintain the metric manifold structure of GR and change the dynamics of the gravitational field 
by effectively giving the graviton a small mass. It is thoroughly assumed that any infra-red (IR) 
modification of gravity must bring in at least one new mass parameter, corresponding to the scale 
at which such a modification becomes effective. However, that is not necessarily the case. In [2] we 
put forward an infra-red modification of GR that does not involve any new adjustable parameter 
but that might require a substantial reassessment of the present understanding of gravity in the 
semi-classical regime. In this paper we review and refine the model sketched in [2] and derive 
some of its cosmological implications. 

We start from usual semi-classical gravity, where matter fields are quantized on a curved 
background manifold, and we modify in the IR the operators of the matter field theory. In 
a cosmological setup, operators corresponding to Fourier modes of physical momentum k are 
corrected by order H 2 /k 2 , H being the Hubble parameter. The coefficient of such a modification 
is fixed by imposing an "Ultra-Strong" version of the equivalence principle, which we are going 
to recall in the following. We argue that such IR-modified operators fail to be supported by a 
metric manifold structure at large scales. Hence, the proposed modification seems to indicate a 
breakdown, for systems of size comparable to the inverse curvature, of the geometrical description 
of space-time as a metric manifold. 

It is known that there is a strict connection between the geometric properties of a manifold 
and the spectrum of the differential operators [3] or the algebra of functions [I] therein defined; 
such abstract characterizations have occasionally been used for generalizing common geometrical 
concepts and the description of spacetime itself [U [5] . However, so far, attempts in this direction 
have always been applied to the UV and intended to modify spacetime at the smallest scales. 
Although perhaps less intuitive, a modification effective in the IR - and at length scales set by 
the curvature - is appealing on both theoretical and observational grounds: 

• The standard paradigm "GR + matter fields" presents some difficulties that are arguably 
UV-insensitive, such as the cosmological constant problem [6] and the black hole information 
paradox [7j. Such difficulties appear already at the semi-classical/low-energy effective level 
and therefore it seems unlikely that any modification at the Plank scale or UV-completion 
of the theory might possibly address them. Questioning the IR-side of the picture seems 
thus a valuable alternative. 

• The apparent acceleration of the Universe is described in GR as something kicking in at 
some late-time epoch in the cosmological evolution and related to an extremely tiny mass 
scale. Such a scale is adjusted ad hoc in the models of "Dark Energy", which, therefore, 
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generically suffer from fine-tuning and "coincidence" problems. The alternative view sug- 
gested here is based on the simple remark that any high redshift object is placed from us at 
distances of order the inverse curvature (i.e. the Hubble radius Hq 1 ). Hence, modifying the 
usual framework at a length scale set by the curvature - rather than fixed a priori - will 
systematically affect any cosmological observation at high redshift, regardless of when such 
an observation takes place and without the need of any external mass parameter. 

Before summarizing the cosmological implications of the proposed deformation let us review 
the main theoretical idea: the Ultra Strong Equivalence Principle [2]. Consider the vacuum 
expectation value of the local energy density of a (massless) field according to standard semi- 
classical gravity; this can be expanded at high momenta as 

<T (M)> bare = Jd 3 k(k + ^^ + ^ + ...^ (1.1) 

= local terms + non local terms. 
Here spatial homogeneity has been assumed for simplicity and the fs are functions of time of 



appropriate dimensions (see eq. (3.17) below, or reference jH], for the explicit expression in the 



case of a massive scalar in a flat FRW Universe). The stress tensor renormalization is a well-defined 



prescription to make sense of the infinities appearing in (1.1 ): the local terms can be subtracted by 
local gravitational counterterms (cosmological constant, Newton constant and higher order local 
operators). The finite non-local pieces left represent the genuine particle/energy content of the 
choosen "vacuum" state. 

It is interesting to note that some of the difficulties afflicting the standard low-energy framework 
for gravity seem to be related to the spacetime dependent subleading terms (the fs in the case 
considered above) in the energy momentum VEV expansion. The black hole information loss 
paradox is related to the presence of the non-local pieces, as they are the ones responsible for 
particle creation and for the evaporation of the black hole. We argue that also the cosmological 
constant problem is related to the presence of those terms, at least if we restrict the analysis to 
the level of non-quantized semi-classical gravity. Such a problem appears in the renormalization 
procedure that we described earlier as the huge amount of fine tuning that is required to match the 



first quartically divergent piece of ( 1.1 ) with its observed value. However, the quartically divergent 
term is not particularly worrisome by itself, essentially because it is a constant. In flat space we 
deal with it, by normal ordering or just by subtracting a constant, albeit infinite, contribution. 
In other words, as brutal as it might appear, "the vacuum does not gravitate" is a well defined 
prescription, as long as the contribution to be neglected is a spacetime-independent scalar. We 



are arguing that the main source of trouble is not the first term in (1.1), but the remaining 
spacetime dependent pieces, as they make the normal ordering prescription meaningless and the 
renormalization procedure necessary in the first place. 



Based on the above reasoning and encouraged by few (more philosophical) arguments explained 



below in Sec. 1.1 , the conjecture at the basis of this paper is that there exists a theory that resembles 



standard semiclassical gravity at small scales but in which, as in flat space, all time- dependent 



pieces of (1.1) (the fs) at any point-event, and for some "vacuum" state, just do not exist. Since 
the terms that we want to get rid of are effectively IR with respect to the first quartically divergent 
part and can be calculated, we should be able to probe this theory by IR-completing standard 
semi-classical gravity. 

It is clearly not guaranteed that a consistent theory with the required properties exists in full 
generality beyond the example considered here; if it does, its formulation will certainly require a 
major theoretical breakthrough. Giving up the metric manifold structure - and dealing only with 
quantum operators - constitutes a significant technical difficulty On the other hand, compared to 
some of the subtleties of standard QFT on curved spacetime [9], the physical picture that we are 
insinuating is somewhat simplified, made trivial. We are conjecturing the possibility of describing, 
locally, the physical phenomena in exactly the same way as we do in flat space, and that the large 
scale-properties of spacetime have no effect whatsoever inside Einstein's free-falling elevator. This 
is the spirit of the Equivalence Principle, that we are upgrading to the level of semi-classical 
gravity and making "Ultra- Strong" [2j. The local resemblance to flat space is pushed to the level 
of the quantized fields, whose energy momentum tensor VEV at any point-event is conjectured to 
amount, for some unambiguously chosen vacuum state, to the constant value ~ J kd 3 k. 

In this paper we try to explore the conjectured theory through a sort of Taylor expansion 
around GR, the zeroth order being a flat FRW Universe with a given scale factor a(t) and a 
spectator free scalar field in it. We isolate the first correction to the scalar field Fourier operators 



that reabsorbs the quadratic divergence / qua d in (1.1 ). In a similar fashion, if we did not know GR, 



we might at least try to explore it by expanding around flat space at a given point in Riemann 
normal coordinates. The proposed modification is first introduced, in Sec. |4| through modified 
commutation relations between the global creators and annihilators of the field, 

(H 2 a 2 \ 
1 - + higher orders . (1.2) 

In the above, n is the comoving momentum and the form of the correction term is dictated by the 
Ultra-Strong Equivalence principle. At this order of approximation, however, there seems to be an- 
other more practical way to see the proposed modification, which we show in Sec. [5j After a simple 
rescaling, one can deal with the usual "manifold-operators" , labeled by the "manifold-momenta" 
n, that satisfy the usual commutation relations, [An, At,] = <^,n'- The proposed modification then 
boils down to introducing a slight mismatch between the "global, manifold" label n and the phys- 
ical (comoving) momentum k that locally defines the infinitesimal translations and the derivatives 



of the local fields: 

S-«(l-^). (1-3) 

Note that the comoving physical momentum is time dependent in this model, although only 
negligibly at sub-Hubble scales. 



The modified dispersion relations (1.3) come from requiring the cancellation of the quadrati- 



cally divergent time-dependent piece / qU ad(£) m the high momentum expansion ( 1.1 ) at some given 
point "x = 0". In Sec. [6] we try to make sense, globally, of such modified dispersion relations 
and explore this spacetime a finite distance away from x = by using the accordingly modified 



translation operator. As much as the physical comoving momentum (1.3) is conserved only ap- 
proximately, it turns out from the analysis of Sec. [6] that the physical distance d between two 
comoving observers grows as the scale factor a(t) only if that distance is much smaller than Hub- 
ble. At large scales, the expansion rate effectively detaches from the local expansion a(t). This 
is the most striking signal of a breakdown of the manifold description. The comoving distance 
A = d/a(t) of two faraway comoving observers, instead of being constant, satisfies 

(H 2 a 2 )' 

A = -^ 3 ^ — - + higher orders. (1.4) 

The effect is clearly negligible at sub-Hubble scales but will affect the observed cosmological 
dynamics at high redshift. 

The cosmological implications/predictions of the present model, obtained by applying the 
proposed IR modification to a spatially flat FRW Universe, are derived and discussed in Sec. [7] 
and can be summarized as follows: 

• Comoving momenta k are conserved only when they are very high compared to Hubble, i.e. 
when k ^> aH. More generally, what is conserved is the quantity 

(H 2 a 2 \ 
1 + -— g — h higher orders J . (1.5) 

• The physical distance d between two comoving observers grows in time proportionally to 
the scale factor a(t) only if the two are well within the Hubble radius. More generally, their 
distances at two different times t' and t are related by 



d{t) _ a{t) 



'l-(t') /„ 2i " 2 t/) u2l 



1 + —- 1 ( H 2 (i)^-j- - H 2 (t') ) + higher orders 



1.6) 



The angular-diameter distance is given, in terms of the redshift, by 
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d A (z) = f^. (1.7) 
1 + z 

In a matter dominated Universe, r(z) is the solution of 

~ltz~ = (1 + ^)3/2 + + hlgher ° rderS - (L8) 

with initial condition r(0) = 0. 

The luminosity distance cIl, at least at this order of approximation, is related to d^ in the 
usual way, di = (1 + z) 2 dA- 



The cosmological distances derived with equation (1.8) in a matter dominated Einstein-de Sitter 
Universe are corrected in the direction of an apparent "acceleration" and show some resemblance 
to those of a ACDM Universe (Figs. [l]and|2j Sec. A more quantitative analysis will follow 

The paper is organized as follows. The general motivations and the main idea are given in 
Sees. | X . X | and |2| which are basically taken from the shorter note [2]. In Sec. [3] we review the 



standard calculation of (1.1) for a free scalar in a flat FRW Universe. In Sec. |4jwe modify the 
Fourier operators in order to get rid of the quadratic divergences there encountered. Modified 
dispersion relations are introduced in Sec. [5] as an essentially equivalent, but more practical, 
way to look at the proposed modification. In Sec. [6] we attempt to make sense, globally, of the 
proposed modification and study its implications on large scales. In Sec. [7] we finally derive 
some cosmological consequences, such as the modified luminosity distance. Potential problems 
and future perspective are finally summarized in Sec. [8] 



1.1 The Equivalence Principle: Gravity as an IR effect. 

The equivalence principle (EP) can be formulated simply as follows: inside a sufficiently small 
free-falling elevator you do not see the (classical) effects of gravity. Amusingly, such a cornerstone 
of modern physics is actually stating what (where) gravity is not, rather than what (where) gravity 
is! Among the many celebrated implications of general relativistic physics, the view that we aim 
to stress here is that EP forces us to consider and describe gravity as an IR phenomenon, whose 
effects are visible only outside the free-falling elevator. How EP turned into a consistent theory 
is well known: gravity is beautifully encoded in GR as the geometry of the physical space-time 
and therefore its effects are automatically suppressed within those systems that are much smaller 
than the inverse curvature. By changing (curving) the large-scale structure of spacetime, GR 
leaves the smallest systems free of classical gravitational effects. Notably, the IR scale where 
non-gravitational physics breaks down is not a parameter of the theory, but is set by the local 
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curvature R. Schematically, in three dimensions, the area of a two-sphere of radius I and volume 
V receives corrections from the flat-space expectation of the type 



A{1) = Anl 2 {l + 0{l 2 R)) = (367r\/ 2 ) 1/3 (l + 0(^ 2/3 ))- (1-9) 

The effects of gravity, originally banned from the free-falling elevator, reappeared, after the 
developments of quantum theory, through what one might call the back door. The fields quantized 
on a curved manifold are sensitive to the global properties of spacetime because their modes are 
defined on the whole of it. As a result, inside the elevator, you will generically experience, and 
possibly detect with your local instrument^] particle creation because of non-local gravitational 
effects. Clearly, the process of quantum particle creation does not contradict EP, which was 
formulated within the framework of classical physics. Nevertheless, it is tempting to try a different 
path than the one historically followed: rather than imposing field quantization on top of a curved 
manifold, here we attempt to upgrade the equivalence principle and extend it to the quantum 
phenomena. Thus, we will consider a stronger version of EP, in which all the effects of gravity 
are definitely forbidden inside the elevator, including the quantum effects that in the standard 
semi-classical treatment lead to particle creation. Effectively, what we want is that for each matter 
field or sector sufficiently decoupled from all other matter fields, there exists a state, the "vacuum" , 
that is experienced as empty of particles by each free-falling observer. 

Following the discussion given in the introduction, a possible concrete strategy would be trying 



to reabsorb in the IR the non-local terms inside (1.1) after the stress energy tensor has been 
renormalized by usual means. However, we will follow a different, more radical, plan. We will 
circumvent the usual procedure of stress-tensor renormalization and just try to cancel the time- 



dependent terms of the bare energy momentum tensor (all terms inside (1.1) except the first). In 
fact, the time-dependent pieces can arguably be reabsorbed with some IR modification, because 
they are effectively IR with respect to the quartically divergent term. Moreover, bypassing the 
usual procedure of stress tensor renormalization might, at the same time, shed some new light on 
the cosmological constant problem, at least at the level of non-quantized semi-classical gravity. If, 
in the IR-completed theory, all the time dependent terms just do not exist, then we do not need 
to renormalize the stress tensor anymore. Of course, we are still left with the leading spacetime- 
constant quartic divergence, but we can live with it and treat it, as we do in flat space, by normal 



1 It is arguable that a small system cannot detect quanta of arbitrarily small energies. While this is generically 
true in practice, there is, to our knowledge, no universal (independent of the details of the matter field dynamics) 
and in principle argument that limits the possibility of particle detection by small detectors. (On the opposite, 
there is a very precise and quantitative sense in which tidal forces in GR become small within small systems.) 
Within a system of size L, the first excited level is typically at energy 1/L above the ground state and cannot be 
excited by photons of much lower energies. However, the structure of the higher energy meta-stable states is often 
very rich, and energy differences among such states can be arbitrarily small. An atom, for example, has extremely 
fine energy transitions that, in principle, allow to "detect" photons of wavelengths much larger than its size. 
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ordering. Finally, it is interesting to reabsorb the quadratic divergence in the IR, rather than with 
a local counterterm, because the required modification, as we are going to show, is of the right 
order of magnitude to give interesting cosmological implications. 

Ultra Strong Equivalence Principle: For each matter field or sector sufficiently 
decoupled from all other matter fields, there exists a state, the "vacuum", for which 
the expectation value of the (bare) energy momentum tensor reads the same as in flat 
space, regardless of the configuration of the gravitational field. 

As a - rather obvious - warning to the reader we mention that, while EP is extremely well 
tested (see e.g. [H]), the proposed ultra strong version is not. Experimental hints might actually 
be pointing against it: the mechanism of quantum particle creation is successfully exploited by 
inflationary models in order to produce the primordial spectrum of cosmological fluctuations. The 
present knowledge of the earliest Universe is, however, too model dependent (and largely based on 
the assumption that such a mechanism is at work) for that to be considered as an experimental 
evidence. 

2 Regions of Space as Quantum Subsystems 

Before describing the model, it is worth examining closely what happens in the standard for- 
mulation of semi-classical gravity, where fields are quantized on a curved background manifold. 
Consider a spacetime with a global time foliation labeled by a time parameter t. The Universe as 
a whole at time t is a three-dimensional manifold M. in the GR description. The matter quantum 
fields are instead described by a quantum state living in a Hilbert space Ti. Thus, if now we 
consider a region of space (at time t) of finite volume V ("this room, now"), that has two comple- 
mentary descriptions [121 [131 E] : it is a submanifold according to GR and a quantum subsystem 
for the quantized fields. The correspondence sub-manifold/sub-system is explicitly realized by the 
set of local operators A(t, x) of the field theory. Those act on the quantum system 7i but have 
labels x living on the three-dimensional manifold M.. As a consequence, we can take integrals of 
scalar local operators over some region of volume V, 



which are still operators acting on the total Hilbert space 7i. More precisely, A(t, V) acts non- 
trivially only on the quantum subsystem Tly corresponding to the region of space that has been 
integrated over, and as the identity on the rest of the system Ti^est- I n this way, the algebra of 
operators A(t, V) determines [T5] the partition H = 7i v ® H^ est of the quantum system/Universe 
and therefore can be taken as a definition of the region of space V. Defining regions of space as 




(2.1) 
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quantum subsystems [TJ] looks like a useless complication in the standard formalism because the 
correspondence sub-system/sub- manifold is always implicitly at work: such a correspondence is 



set once and for all by (2.1 ), given the set of local operators A(t, x). 

When we say that the metric-manifold description might break-down, we mean, more precisely, 
that regions of space, still perfectly defined as quantum subsystems, may not have the nice property 



that the corresponding operators integrate as in (2.1). We argue therefore that (2.1) is valid only 



in the limit of zero curvature and that operators corresponding to extended regions of space can 



be written as (2.1) only up to order 0(RV 2 ^ 3 ), where R is some curvature scalar. Here is the 



mnemonic thumb-rule that will guide us in our further developments: 

A(t, V) ~ / d 3 aV-0 (3) A{t,x) [l + 0(RV 2/3 j] . (2.2) 
Jv 



We are deliberately mimicking the general type of corrections (1.9) that non-extensive geometrical 
quantities undergo in the transition from flat space to curved space. The idea is to extend this 
type of behavior also to extensive quantities, that in the standard description are just proportional 



to the volume. The implications of (2.2) are quite striking. Consider, for instance, a perfectly 



homogeneous Universe. By definition, each comoving observer measures in its surrounding the 



same energy density p. Eq. (2.2) implies that in that Universe, if one starts considering regions 
of space of Hubble size, the total energy inside that region will drastically differ from the three- 
dimensional integral of the local densities measured by the observers living therein. In a metric 
manifold, that would be a dramatically non-local effect. 



3 The Zeroth Order: a Flat FRW Universe 

In this section we briefly review the quantization of a free scalar field in a FRW spatially flat 
metric manifold. In the spirit of this paper, what we are considering here is therefore the zeroth 
order approximation of the theory that we are trying to explore. From the action 

S= \ I d 4 x^ (<90 2 - m 2 2 ) (3.1) 



2 _ 

and the metric element 

ds 2 = dt 2 - a 2 {t)(dx 2 ) = a 2 (r)(dr 2 - dx 2 ) (3.2) 
we derive the equation of motion for the field, 

d 2 

0(t, x) + 3H<j)(t, x) x) + m 2 (j)(t, x) = 0, (3.3) 



a{t) 



and its energy momentum tensor, T^ u — 1 

id comoving s; 

T °(t,£) = H(t,x 



— S J^ , from which we get the expression of the energy 
density of at time t and comoving spatial coordinate coordinate x: 



2 + ^0) 2 + m 2 2 



(3.4) 



A dot means derivation with respect to the proper time t while a prime denotes derivation with 
respect to conformal time r; H = a/ a is the Hubble parameter. Upon quantization in the 
Heisenberg picture, the field can be expanded in creators and annihilators: 



(f){t,x) 



(2tt) 3 / 2 



drf 



(p)t 



(3.5) 



For reasons that will be clear in the following, we have called n (and not k) the comoving momen- 
tum label in this FRW space; n is a conserved quantity, related to the proper physical momentum 
pbyp = n/a(t). The suffix (0) stands for "zeroth order" since the creation operators just intro- 
duced are intended to be the zeroth order (-metric manifold) approximation of the operators of 
the IR-completed theory. 

From the action (3.1) one derives the conjugate momentum 

_ 5S 

5cf){t,xy 

and from the canonical commutation relations 

[4>(x), tt(x')] = i5 3 (x - x'), 
one derives the canonical commutators for the zeroth order- global operators 

[Af,A^) = 8\n-n>). 



(3.6) 



(3.7) 



(3.8) 



It is customary to choose A% as the operator that always annihilate the vacuum. The chosen 
vacuum state is therefore implicitely characterized by the choice of the mode functions ip n (t), that, 
by (3.5) and (3.3), satisfy 



where 



ljj n + 3Hljj n + Lofyn = 0, 



+ m z 



(3.9) 



(3.10) 



There is an extensive literature about the correct vacuum choice, i.e., the appropriate solutions 
of (3.9). The classical argument in favor of the adiabatic vacuum [HI dSl IB] is its resemblance 
to the flat space behavior for those modes that are deep inside the Hubble scale. Within the 
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standard framework, the adiabatic vacuum is as close as we can get to implement the general 
ideas described in the previous sections and, therefore, will be our starting point. The adiabatic 
vacuum corresponds to the choice 



1 



exp 



W n dt' 



(3-11) 



y/2a 3 W n 

where the function W n (t) can be approximated at arbitrary adiabatic order with a formal WKB 
expansion. For the next to leading order in the momentum n expansion we need no more than 
the next to leading adiabatic order [TBI |8] : 



W n (t)=uj n (t)(l + e 2 ) 1 / 2 , 
where the function e 2 can be derived from a WKB expansion and amounts to 



d 



f2 



(a 3 uj n ) 3 / 2 dt \_uj n dt 



1 d 



2 

a I a 



n 2 \a 



+ BT' 



1 + 



m 



rr 



(3.12) 



(3.13) 



The vacuum expectation value of the energy density at a given point is calculated straightfor- 
wardly: 



(o|t °M)|o> 



2(2vr) 
1 

47T 2 



d 3 n d 3 n' 



+ rn 



n ■ n 



n 2 dn(\ip n \ 2 + ujl\ip n \ 2 ), 



[4°US (3-14) 
(3.15) 



where commutation relations (3.8) have been used in the second line. From (3.11) and then (3.12) 
we get 



(\i} n \ 2 + UJ 2 n \*lj r , 



2a 3 W, 
1 



ul + W 2 + - (In a 3 W n y 2 



(3.16) 



a 3 u n 



ul + - 3H + 



+ 0(n 



-2\ 



Note that e 2 drops from the expression of the energy momentum tensor VEV if we do not consider 
beyond the quadratic divergence. In other words, all we need of W n is its expression at leading 
adiabatic order, W n ~ u n . We finally get 

H 2 a 



T °(t,f)|0) 



An 2 a 3 



to be opposed to the flat space result 

(0|T °(t,x)|0) flat 



n dn 



1 



u n + 



2n 



+ 0{n 



47T 2 a 3 
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n 2 dn uj n . 



(3.17) 



(3.18) 



4 Getting rid of Quadratic Divergences in a toy-model 
Universe 



In order to enforce the Ultra Strong Equivalence Principle we want to get rid of the terms in the 
local energy momentum VEV that we normally would not find in flat space. In this paper we deal 
only with the time dependent quadratic divergence, the second term inside the square brackets in 



(3.17). 



First of all, we focus our attention on a given comoving observer /trajectory at any proper 
time t and, say, at point x — 0. For brevity, when we write x ~ 0, we refer to a region of space 
around that trajectory small enough that the manifold description holds and we can define spatial 
gradients and local commutation relations between neighboring points. 

We want to strictly preserve the local dynamics of the scalar field. Our starting points will be, 
therefore, the field equations (3.3) in the Heisenberg picture for the field operator in x = 0, 



<j>(t, x « 0) + 3H<j>(t, x « 0) + 



m 



(f)(t,x « 0) = 



(4.1) 



and the Hamiltonian density T~C(t, x ~ 0) at point 0, which, again, is an explicitely time dependent 
operator: 

di<p 2 (t,x « 0) 



T U = H(t, x^0) 



\t,x^0) + 



+ m I (t)\t,x « 0) 



(4.2) 



While (3.3) and (3.4) are both elegantly derivable from the same action, for the moment we are 
forced to consider (4.1) and (4.2) as two basic - and in principle independent - ingredients of this 



novel framework. Note moreover that, in a sort of pre-geometric fashion, we can take the above 
equations as the implicit definition of the local expansion a(t). By homogeneity, we will require in 



Sec. 6 that the translated fields at some other point in space still obey (4.1) with the same scale 



factor a(t). Commutation relations between local quantities are strictly preserved: 

[(f)(x), ti(x « x)] = i5 3 (x — a?), 



(4.3) 



(the momentum conjugate to <fi can still be defined completely locally in terms of the Lagrangian 
density in x = 0: n(x w 0) = dC(x w Q)/d<fi(x w 0)). 

As before, we formally expand the field in Fourier modes: 



<f>(t, x w 0) 



(2tt) 3 / 2 



(4.4) 



where the suffix (1) means that we now want to consider the modes corrected at first order in 
the IR expansion. So far the field is only defined at x = 0; the x labels at the exponent in the 
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equation above are there just to remind that, if we take space derivatives of the such a local field, 
a factor in drops. The crucial observation now is that the commutator between A and A* is 



proportional to the total volume and therefore will receive the postulated corrections (2.2). By 



looking at (3.14) and (3.17), it is immediate to note that the quadratic divergence is reabsorbed 



if we impose the commutation relation 



[^UrH^-n'lll-^]. (4.5) 



The correction is in fact of the expected form (2.2). This is perhaps more evident in the approach 



followed in [2] where a compact toy-Universe of total linear size L much smaller than the Hubble 
length was choosen to begin with. In that framework, Fourier labels assume discrete values and 
the correction H 2 a 2 /n 2 is clearly of order H 2 L 2 or smaller. Although we previously referred to 
modifications at scales of the inverse curvature, we find that, more precisely, it is the extrinsic 
curvature which appears to regulate the correct modification. Equivalently, we can just say that 
the annihilators are modified with respect to GR by 



The corrected Fourier modes are effectively time dependent. At least at this order of ap- 
proximation, however, the time dependence is just a multiplicative factor that leaves the vacuum 



untouched. Note also that, by applying the postulated local equation (4.1), the time dependence 
of An^ introduces additional spurious terms inside (3.9). Those, however, are of higher order in 



the n expansion and therefore do not affect the present result. Finally, inverse Fourier transforms 



(such as inverting (4.4) by integrating over x) are not allowed. Following the reasoning of section 
[2J we argue that whenever we take integrals J d 3 x of operators over a volume V we are making 
an error of order 0(H 2 V 2 ^ 3 ). 

So far we concentrated only on local operators at some given point x = and made connection 
with the global Fourier modes. A crucial ingredient to understand the large-scales properties of 
this deformed spacetime is the momentum operator constructed with the modified Fourier modes, 

pW= [ d 3 nnA^A^. (4.7) 



In fact, we are not provided with a metric element that we can integrate over, but we can exploit the 
assumed homogeneity of this spacetime and define a translation operator by simply exponentiating 
(4.7). By applying the translation operator, in section [6] we are able to define the local operators 



a finite distance away from x — 0. 
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5 Refining the Model 



The form of the momentum operator suggests another way to look at the proposed deformation. 
We can write (4.7) in terms of the usual GR operators A^\ satisfying the usual commutation 



relations, and modify the momentum labels instead: 

p (1) = |rf 3 nfc4 0)t 4 0) , (5.i) 

where, at this order of approximation, 

* = a(i__j. (5.2) 

We are trading modified Fourier operators for modified dispersion relations. In other words, in 
this equivalent, but more practical, formulation, we are postulating a slight mismatch between 
the "metric-manifold" labels n and the physical momenta k that locally define the infinitesimal 
translations and the derivatives of the local fields. Since n is the usual comoving "manifold" label, 
we assume that it is time-independent. Moreover, since we will use only the usual "manifold" 
creators operators from now on, we drop the the suffix (0) to simplify the notation. The 
recipe can finally be summarized as follows: 

• There are two different (comoving) Fourier coordinates: the "manifold" -coordinate n and 
the locally defined physical momentum k. The integration measure in Fourier space is flat 
in the n-coordinates. Moreover, during time evolution, ft is conserved, k is not. The relation 
between the two is 

(H 2 a 2 \ 
1 ^ 2~ + higher order j . (5.3) 

• The local field in x ~ is expanded in Fourier modes as follows: 



1 



cf>(t,x « 0) = j^j- 2 J dn A ^ k {t)e ik '*An + r k {t)e- ik - s A\\ ee J dn^e** 3 (5.4) 

and the A n satisfy the usual commutation relations: 

[A H ,Al} = 6 3 (n-n'). (5.5) 



Again, the x coordinate in (5.4) extends only as far as is needed to take the space derivatives 
of the field in x = 0. When a derivative is taken, a factor of k, instead of n, drops. 
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• The local equation for <p(t,x ~ 0), (4.1), fixes the equations for the functions ip n : 

4> k (t)+3HiP k (t)+co 2 k ip k (t) = 0, 



(5.6) 



where 



k 2 



+ m 2 



(5.7) 



and k is given in ( 5.3[ ). The solution of the above equations are choosen to correspond to 
the adiabatic vacuum (in the physical momentum k), i.e. 



exp —i 



W k dt' 



(5- 



yWV k 

The function W k can be approximated at arbitrary adiabatic order as before. The only 
difference is that u k should substitute uo n in (3.12) and (3.13). 

The energy density in x ~ is expressed in terms of the local field 4>(t, x ~ 0) in the usual 
way (equation (4.2)), 



T ° = H{t, x « 0) 



<fi 2 (t,x « 0) + 



d i( j) 2 (t } x « 0) 



+ m 2 (f) 2 (t } x « 0) 



(5.9) 



The above equations can be used to calculate the vacuum expectation value of Tq. By direct 



substitution of (5.4) into (5.9) we get 



(0|I*(t,a?«0)|0) 



1 



2(2tt) 3 

J /"OO 

4^2 



d 3 n d 3 n' 
n 2 dn (|^fc| 2 + ^|^ fc | 2 



(5.10) 



By substitution of (5.8) we then get 



(0|T °(t,x«0)|0) 



An 2 a 3 



4iT 2 a 3 



n dn 



H 2 a /n/ 



n 



''dn [u n + 0(n- 3 )] 



(5.11) 
(5.12) 



Note that the quadratic divergences are here reabsorbed just be re-expressing (0| Tq |0) in terms 
of the appropriate "flat measure" time-independent Fourier coordinates n. 
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6 A look at the global picture 



So far, we have defined local operators only in x ~ and made connection with the global ones. 
However, we want to be able to move in this space arbitrarily far in any direction and define local 
operators also elsewhere. We are provided with a momentum operator, 



where, again, 



Pi 



k(n) 



d 3 n k(n) AlA H 



ft 1 



H 2 a 2 
2n 2 



(6.1) 



(6.2) 



Since we are assuming homogeneity, we obtain a global translation operator by exponentiation, 



T i (\) = e~^. (6.3) 
Thus, we can define local operators at a finite distance from x ~ by the action of Tf. 



<f>(t,*)=Ti(\) <f>(t,0)Tr\X) 



(2vr) 3 / 2 



d 



3 i -iXnJl-S^) 



(6.4) 



where * stands for "at distance A if one keeps going in the i direction" and A is made up by 
all the infinitesimal steps that have been taken in order to go from one point to another and is 



therefore the (comoving) proper distance. The momentum operator is given in (6.1)-(6.2) only 



approximately. Therefore, also (6.4) should be trusted only as far as d = a\ is smaller than the 
Hubble radius. 

In the present case, because of homogeneity, translations in different directions commute with 
each other, so we can label the local operators with three parameters and describe, approximately, 
a finite patch of Universe around x ~ 0. Explicitely, 



<f>(t,X) 



d 3 n 



(6.5) 



(2vr) 3 / 2 

Although the above equation has a familiar form, the labels A* should really be thought here as 
the parameters of the translation group. We will soon show that, if we use them as integration 
variables at t — const., as we would do on a standard FRW flat space, we make an error of order 
V 2 / 3 H 2 . 

It is interesting to calculate, at some given time t, the canonical commutator between 7r(0) = 
a 3 (p(0) and the field itself at a distance. By using the standard commutation relations ( 5.5[ ) we 
obtain 



[7T(0), 0(A)] 



(27T) 



d 3 n 
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ik-X 



(6.6) 



By using the Wronskian condition a 3 (iplifjk — tfk'&k) = immediately derivable from (5.8), and 
(6.2) we get, up to higher orders in Ha\, 



[tt(O),0(A)] 



d ne 



-i A-n( 1- 



(2n) 



d I n e" 



i X-n 



H 2 a 2 
2n 2 



1 + iX ■ ft 



(6.7) 



where we have expanded the exponential in the last expression. The first term gives the usual 
delta function. The last term gives 



iH 2 a 2 d 



d 3 



n 



-i a X-n 



a=l 



i H 2 a 2 

87T A 



2{2nfda 

In this way we obtain the first correction to the canonical commutator, 

[tt(O),0(A)] = -z^ 3 (A) + 1 Hl(r 



57T 



A 



(6.8) 



(6.9) 



The above relation finally gives a more precise sense to (4.3). At face value, (6.9) is a signal 
of a non-locality that becomes more and more important at large scales. But again, that is the 
interpretation that we would have of (6.9) on a metric manifold. In section [2] we argued that 



regions of space at a given time can still be associated with a well defined algebra of "averaged" 
operators of the matter theory, but that those operators can only approximately be expressed as 
integrals over a metric manifold. The expression in (6.9) is in fact a three-dimensional density 
and the correction term is just signaling the mistake that we make if we naively take the integral 



ci 3 A'[7r(O),0(A')] 



Xa<^H~ 



1 + i#VA 2 + 0{H 2 V 2 ' 3 ) 2 



(6.10) 



and use A as a flat-measure integration variable. We nicely recover the expected pattern (2.2). 



It is equally interesting to calculate [0(0), 7r(A)]. First, we note that there is a potential 
ing tl 

tt(A) 



ambiguity in defining the time derivative of a displaced operator. By deriving (6.5) we get 



(2vr) 3 / 2 

^3 



d A n 



-ik-X 



+ 



d 3 n (k ■ A)' 



—ik-X 



At 



(6.11) 



(2tt) 3 / 2 

The second line in the above equation is there because k is time dependent. However, if we just 
applied the translation to vr(0), instead of deriving 0(A), those terms would not be there. Therefore, 
for consistency, we need to make them ineffective at the required order of approximation. In other 
words, we have to impose that [0(O),7r(A)] = — [7r(0), 0(A)]. Because of the second line of (6.11), 
the commutator between 0(0) and vr(A) gives 

^3 



[0(0), tt(A)] = -[tt(O), 0(A)] -2i 



(27T)* 



d i ne~ ik ' x \ip n \ z (k ■ A)' 



(6.12) 
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the last term being the spurious contribution. In order to get rid of it, we have to make the co- 
moving physical distance A also time dependent. This effectively means that, after an infinitesimal 
time step dt, we have to reconsider the field translated, from x ~ 0, not by the same comoving dis- 
tance A, but by a slightly different amount. At high momenta/small distances, since \ipk\ 2 ~ I/ 1 "-, 



the integral in the second term of (6.12) reads 



d 3 ne 



-in-X 
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A • n 1 



H 2 a 2 



A • n 



(H 2 a 2 
2n 2 



(6.13) 



We make the ansatz A = /3A (H a )'A, (3 being a number to be determined. We get 



d 3 ne- iU \^ n \ 2 {k-\) ~ i{H 2 a 2 )^ [ d 3 n ( ^ 



1 



—ia n-X 



(6.14) 



a=l 



The last integral can be regularized by setting n 



n 2n 3 

-•! , n~ 3+e and taking the e — > limit only after 
deriving with respect to a. The result is null for (3 = 1/4, which fixes the time dependence of A: 

< ^(H 2 a 2 ) , . 

^ = ^ — - A — ■ (6-15) 



The above formula is one of the main results of this paper and can be rephrased as follows. 
In this homogeneous Universe, in the vicinity of each point, there is a local expansion given by 
the same scale factor a(t). In essence, a(t) rules the evolution equations for the local operators 
and determines the rate at which the distance between two comoving observers grows, as far as 
such a distance is extremely small compared to Hubble. In that limit, the comoving distance A is 
conserved. However, if you pick up a pair of comoving observers further apart, their distance grows 
slightly differently, which means that the expansion rate on large scales is effectively dependent 



on the distance. The modified expansion rate is given, at first order, by (6.15). This effect, which 



is the most striking signal of the breakdown of the manifold description, is clearly negligible well 
within the Hubble scale. 



7 Cosmology 



Cosmology is where the predicted departures from GR are more relevant, since every high red-shift 



object is located from us at distances comparable to the inverse curvature. Eq. (6.15) can be 
integrate straightforwardly: 



X- 2 (t) - \- 2 (t>) = - l -[H 2 (t)a 2 {t) - H 2 (t')a 2 (t')]. 
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(7.1) 



Equivalently, the proper distance d 
t as 

d(t) _ a(t) 
Kf) ~ a(f) 



aX between two comoving observers evolves from time t' to 



1 - 



d 2 (t') 



H\t) 



a 2 (t) 



H 2 (t') 



1-1/2 



(7.2) 



2 V" V a?{t') 

According to (7.1 ) and ( |7.2 ) the distance between two comoving observers grows slightly less than 
in a FRW Universe with scale factor a(t). 

We now want to calculate the "corrected" trajectory r(r) for a light ray, where r is its comoving 
distance and r the conformal time. A light ray passing through x ~ satisfies dr = dr. However, 
on top of the usual contribution, the trajectory of a light ray receives a correction from the modified 



global expansion (6.15): 



dr 
dr 



1 + r' 



(7.3) 



Now we specialize to a matter dominated Universe for simplicity. Therefore, we set Ha = 2/r 
and we get 



dr 
dr 



(7.4) 



Normally, at the RHS we would simply have 1 and, again, the second term should be considered 
just as the first order correction. We couldn't find an analytic solution to the above equation, and 
therefore we will solve it numerically. With respect to an Einstein-De Sitter Universe, a light ray 
propagating forward in time will go a shorter way. A light ray that we are detecting now, on the 
other hand, is actually coming from farer. For a light ray propagating back in the past the above 
equation becomes, in terms of the redshift, 



d(H r) 
dz 



1 



(1 + zfl 2 



+ 



(H or y 



(7.5) 



The above should be solved with initial condition r(0) = and will give some r(z) to be used in 
what follows. Note that, as in the normal framework, the Hubble parameter H is just a total 
normalization factor for cosmological distances. 

Light rays going backward in time in any direction describe, at some redshift z, the surface of 
a two-sphere of area 4Tra(z) 2 r(z) 2 . Therefore, the angular diameter distance is r(z) multiplied by 
a(z): 



d A (z) 



1 



-r(z). 



(7.6) 



l + z 

The redshift of the frequences, on the other hand, is basically unmodified in this framework, 
1 + z = 1/a. It is true that the comoving momentum k (6.2) is not strictly conserved, but the 
effect is negligible for momenta much bigger than H . 

Obtaining the formula for the luminosity distance looks, in the present framework, slightly 
less straightforward. We discuss its derivation in the Appendix where we argue that, at least of 



It 



(a) Present Framework 



(b) Einstein-de Sitter 



Figure 1: In both figures the solid curves show the angular diameter distance d,A for ACDM model 
with fl m = 0.3 (upper curve) and for Einstein-de Sitter (Q m = 1, lower curve) where the Hubble 
constant H has been normalized to 1. In the left panel we compare that behaviour with what 
we find in the present model. In order to obtain curves similar to ACDM we have to lower the 
Hubble constant Ho which is the only parameter left free in this model. From top to bottom we 
show the results for Ho = 0.78,0.82,0.86. The runaway behaviour after z ~ 2 is probably due to 
the fact that we are getting very close to the Hubble radius at that redshift (we find, at z = 2, 
d,A — (2/3)/f _1 ) and therefore higher order IR corrections to our formulas are needed. The right 
panel shows the effect of lowering the Hubble constant to those same values in the unmodified 
matter dominated Einstein-de Sitter Universe. 

this order of approximation, angular diameter and luminosity distances are related by the usual 
formula 

d L = (l + z) 2 d A . (7.7) 

We plot the obtained luminosity distance as a function of the redshift in Fig. [2] for different values 
of the Hubble parameter. 

Our correction makes cosmological distances larger than in a usual Einstein-de Sitter Universe. 
This goes in the right direction of a ACDM Universe, at least qualitatively (Figs. [T] and [2j left 
panel). However, this effect kicks in only at a relatively high redshift. By expanding r(z) as 
r(z) = riz+r 2 z 2 + . . . it is easy to see that the proposed effect only modifies the curve at forth order 
in the small-z expansion, while dark energy models usually give already contributions at second 
order. Qualitatively, in order to raise the overall angular diameter distance to the plateau reached 
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(a) Present Framework 



(b) Einstein-de Sitter 



Figure 2: Same as Fig. [I] but for the luminosity distance = (1 + z) 2 d,A- 

by ACDM, we have to lower the Hubble constant H , so that cIa(z) starts out with a slightly higher 
derivative. This makes our curves systematically higher than ACDM at small redshift. Note, on 
the other hand, that lowering H in the absence of this effect produces curves of a completely 
different shape (Figs. [l]and|2j right panel). Finally, the runaway behaviour of d,A after about z = 2 
(Fig. [TJ left panel) is probably due to the fact that the distance covered by the light ray becomes 
dangerously comparable to Hubble at that redshift (indeed, (1a{z = 2) ~ (2/3)if~ 1 (2; = 2)) and 
therefore higher order corrections to the present model are needed. 

While the low-redshift discrepancies with the successful concordance model suggest that this 
mechanism might not be able, by itself, to give full account for the apparent acceleration of the 
Universe (for a more quantitative analysis see PH]), it is worth noting that the case for an Einstein- 
de Sitter Universe with a lower Hubble constant has occasionally been made in the literature. 
Interestingly, the authors of p2] claim that a model with Q m = 1 and Hq = 46 km/s/Mpc fits 
cosmological observations, with the exception of supernovae, better than the concordance model. 

8 Discussion 

Within the established low-energy effective framework "GR + matter fields" 

S = J {M&R + £ ma ttcr) (8.1) 
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there are models that can fit our cosmological data with impressive accuracy. However, those 
explanations do not come without a cost: we need to assume two epochs of accelerating expan- 
sion (inflation and "dark energy"), characterized by widely separated mass scales, that require 
appropriate negative pressure components and, expecially for dark energy, a tremendous amount 
of fine tuning. With the possible exception of inflation, these difficulties look oblivious to the 
possible UV-completions of (8.1). The candidate theories of quantum gravity, as long as they 
reproduce (8.1) at low-energy, are not expected to essentially simplify the picture, although they 
proved able to suggest new mechanisms, give insights on the cosmological singularity, motivate 
the introduction of new degrees of freedom/sectors needed for the acceleration etc. . . 

In this paper we questioned the IR-side of the paradigm "GR + matter fields" and attempted 
to look at (8.1) as a small distance approximation. The proposed modification is radical because 
it challenges the usual geometrical description of spacetime on the largest scale, but compelling 
because it does not contain any adjustable parameter: in a cosmological set-up deviations from 
GR are calculable and become relevant, at any epoch, at the Hubble scale. Abandoning the well 
understood and reliable framework (8.1 ) is adventurous to say the least; but we tried to do that by 
following a physical principle that seems to address, at least in part, some theoretical difficulties 
such as the cosmological constant problem and the black hole information paradox and, at the 
same time, to give observable cosmological predictions at high redshift. 

The most relevant cosmological effect is a departure at large scales from the local expansion 
rate (e.g. eq. (1.6)). The general picture that this model suggests is that the present acceleration 
of the Universe may not be due to any additional "dark energy" component but to a systematic 
breakdown of the geometrical description of spacetime at length scales comparable to Hubble. 
The proposed effect does not kick in at some given time "by coincidence" during the cosmological 
evolution, but is permanently present at any cosmological epoch whenever a comoving observer 
looks at very far away objects. 

Exploring new directions often opens up potential problems and questions which have been 
answered here only in part and surely need future work and further double checks. We isolated 
the IR corrections that apply to a free scalar field in a spatially flat FRW Universe. It is not 
clear if those corrections have some sort of universality i.e. if they automatically apply to, or are 
inherited by, other types of fields. For instance, in a cosmological set-up, applying the Ultra-Strong 
Equivalence Principle to the electromagnetic field is just trivial, because of conformal invariance. 
For photons, the time dependent terms that we got rid of in the high momentum expansion (1.1) 
are absent from the very beginning. If the proposed model is really a modification of the underlying 
geometry, such a modification should be universal, and this, at present, is not guaranteed. 

It would be interesting to consider the implications of this model in the very early Universe and 
close to the singularity, where IR corrections of order Hubble might become important. We note 
also that the usual cosmological kinematic of the modes ( "exiting" during acceleration, "entering" 
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during deceleration) is modified in the present framework, as the physical comoving momentum 
is not conserved on large scales (eq. 1.5). It would be interesting to see what are the implications 
of this effect in the early Universe, although higher orders in the IR expansion might be needed 
in order to have a more complete picture. 

Whether or not the proposed modification applies to reality, our findings are of some interest 
and can be rephrased as follows: how much should we modify GR in order to "renormalize" the 
stress tensor (1.1 ) from the IR side - with an IR modification - instead than with local gravitational 
counterterms? How different from GR would the conjectured theory be? We have shown that the 
required modification is indeed extremely tiny. For a free scalar in a cosmological set-up, Fourier 
modes of physical momentum p get modified by an amount of order H 2 /p 2 . On cosmological scales, 
such a modification goes in the direction of an effective acceleration but is still too weak to give 
full account of supernovae observation [10] . It will also be interesting to modify, along these same 
lines, a black hole or a stellar solution, also to check the compatibility with constraints coming 
from the solar system and strongly self gravitating objects. At the level of order of magnitude 
estimate, however, we do not expect relevant corrections. The proposed modification affects those 
phenomena whose relevant length scales are of order the inverse curvature. Atomic spectra and 
and nuclear physics should be untouched, even in the presence of strong gravitational fields. The 
reason is that the inverse (Ricci) curvature is roughly of order M^/p, p being the energy density, 
and this typically gives a huge number in atomic size units. For what concerns the solar system, 
its entire dynamics is well within the average curvature radius. If we try to compare the typical 
solar system distances to, say, the Weyl curvature, those are always extremely small. An order 
of magnitude estimate gives (Weyl curvature) _1 / 2 ~ r(r/3Km) 1 / 2 , r being the distance from the 
sun. This means that, within the solar system, and as opposed to high-redshift cosmology, objects 
are always at distances much smaller than the inverse curvature. 

We finally note that some recent observational inputs in cosmology may look as an encour- 
agement for IR modifications of gravity. Beside the acceleration of the Universe, which is already 
puzzling by itself, some emerging tensions between standard ACDM cosmology and large-scale 
observations might be indicating that there is something going wrong with the standard GR de- 
scription of the largest scales. Interestingly, other IR modifications of GR have already been 
tested [El [19] against those discrepancies, which include, among other, the lack of large-scale 
CMB correlations (20], the anomalously large integrated Sachs- Wolfe cross- correlation [21] and 
the abnormally large bulk flows on very large scales 
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A The luminosity distance 

The luminosity distance accounts for the flux of energy per unit time arriving from the source. 
In standard cosmology it is known that there are two in principle-independent dampening cos- 
mological effects for the flux. The usual analysis nicely relies on the wave-particle duality: First, 
each photon has its energy red-shifted by a factor 1/(1 + z). Second, the number of them per 
unit time decreases, from emission to absorption, by the same factor. As already discussed, the 
present framework does not change the first effect: comoving momenta much higher than Hubble 
are conserved and therefore physical momenta redshift in the usual way. Here we argue that 
also the second dampening effect is proportional to 1/(1 + z), at least at the present order of 
approximation. In principle, this is not so obvious. We might in fact reason as follows. 

Consider a bunch of photons leaving the source at redshift z and propagating towards us. 
Say that, at a given redshift z, the comoving distance between the first photon and the last is 



A(z) = ra vst (z) — ri ast (z), where both rfi rst (-2) and ri ast (,2) are solutions of (7.5). In the limit in 



which those two solutions are very close to each other it is easy to estimate A(^). By writing 



equation (7.4) for r first (r) we get 

2 



r( ast + A' = 1 - ^(rf ast + 3r? ast A + 3r last A 2 + A 3 ). (A.l) 
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We can now use the equation for ri ast and drop higher orders in A. This gives an estimate for a 
small spread A (z) around a given solution r(r), 

A' r 2 (r) /A . 

By reconverting in redshift this gives 

A(z) = e ^J?^V). (A.3) 

Numerical integration shows that the above is indeed a very accurate estimate. This way of 
reasoning would suggest that, on top of the usual spreading effect due to the expansion, there is a 



sort of "refocussing" effect, due to the fact that two infinitely close solutions of (7.5) tend to get 
closer and closer to each other as they evolve forward in time. For faraway objects that is a very 
relevant brightening effect. 
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What is not convincing about the above derivation is that it heavily relies on the estimate of a 
small spread A even when that is calculated very far away from us. We know that the translation 



operator (6.3) that we are using is very accurate only in shifting objects relatively close to the 
center x ~ 0. The error in estimating a small shift A(z) at high redshift can be comparatively 
very large. 

A more convincing way to proceed is simply to note that two parallel light rays, in the limit 
when they are very close to each other, should feel only the local expansion a(t). At each time 
step, by homogeneity, we can put one of the two at x ~ and calculate the position of the other. 
This amounts to re-set, at each time step, ri ast = in eq. (A.l) and gives, for the comoving spread 



A, an equation identical to (6.15). The effect is negligible each time the spread is much smaller 



than Hubble. We therefore obtain, at least at this order of approximation, the known relation 



di 



[l + z) 2 d A . 



(A.4) 
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